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Abstract 

Using an approach developed in physics, we propose a new framework for the study of cellular networks. The key idea of 
the physical network model we propose is to replace the discrete base stations (BS) entities by a continuum of transmitters which 
are spatially distributed in the network. This allows us to establish a closed form formula of the other-cell downlink interference 
factor /, as a function of the location of the mobile. We define here / as the ratio of outer cell received power (i.e. the power 
received from other cells) to the inner cell received power. This physical model allows calculating the influence of interference 
on any mobile in a cell, whatever its position. Results obtained with that closed-form formula are close to the ones obtained by 
simulations using a traditional hexagonal network model. Since the physical model allows to establish a closed form formula 
of the interference factor, it allows to do analytical studies of wireless networks such as outage probability, quality of service, 
capacity. 

I. Introduction 

The estimation of cellular networks capacity is one of the key points before deployment and mainly depends on the 
characterization of interference. As downlink is often the limited link w.r.t. capacity, we focus on this direction throughout this 
paper, although the proposed framework can easily be extended to the uplink. In CDMA and OFDMA systems, an important 
parameter for this characterization is the other-cell interference factor / (OCIF). The precise knowledge of the OCIF allows 
the derivation of outage probabilities, capacity evaluation and then, the definition of Call Admission Control mechanisms. 

On the downlink, [1] [2] aimed at computing an average OCIF over the cell by numerical integration in hexagonal networks. 
In [5], Gilhousen et al. provide Monte Carlo simulations and obtain an histogram of /. In [4], other-cell interference is given 
as a function of the distance to the BS thanks to Monte-Carlo simulations. Chan and Hanly [6] precisely approximate the 
distribution of the other-cell interference. They however provide formulas that are difficult to handle in practice. 

In contrast to previous works in the field, the modelling key of our approach is to consider the discrete BS entities of a 
cellular network as a continuum. Recently, the authors of [11] described a network in terms of macroscopic quantities such as 
the node density. The same idea is used in [12] for ad hoc networks. They however assumed a very high density of nodes in 
both papers. 

This paper extends the framework proposed in [7] and [9] that provides a simple closed-form formula for / on the downlink 
as a function of the distance to the BS, the path-loss exponent, the distance between two BS, and the network size. We show 
that a network can be considered as a continuum of BS, by using an analogue approach as in physics. 

The rest of the paper is organized as follows. In Section [II] we first introduce the interference factor /. It characterizes the 
relative influence of the interference generated by transmissions originating from the other cells of the network. In Section [Hi] 
following an analogue approach with the gravitation's case, we show that a transmitting BS can be replaced by a density of 



BS distributed all over the cell. In Section IV we generalize that approach to a network. Considering the base stations of a 
network as a continuum, we establish a closed form formula of the interference factor /. . By comparison with simulations, 
we establish that the physical model and the traditional hexagonal model are two simplifications of the reality. None is a priori 
better than the other but the latter is widely used, especially for dimensioning purposes. That is the reason why a comparison 



is useful. We show that the physical model is accurate even when the density of BS is very low (see section IV-B i In Section 
V we conclude. 

II. Interference model and notations 

We consider a cellular radio system and we focus on the downlink. BS have omni-directional antennas, so that a BS covers 
a single cell. If a mobile u is attached to a station b (or serving BS). The following power quantities are considered: 

• Pb >u is the transmitted power from station b towards mobile u (for user's traffic); 

• Pi, — Pcch + ^uPb.u is the total power transmitted by station b. In CDMA systems, P cc h represents the amount of power 
used to support broadcast and common control channels. 

• pt.u is the power received at mobile u from station b; we can write pb u = u ; gi, u designates the pathloss between 
station b and mobile u 

• Sf,,u — Pb,u gb.u is the useful power received at mobile u from station b (for traffic data); 

The total amount of power experienced by a mobile station u in a cellular system can be split up into several terms: useful 
signal (S u ), interference and noise (Noise). It is common to split the system power into two terms: I u = Ii n t,u + Iext,u, where 



Iint,u is the internal (or own-cell) received power and I ex t tU is the external (or other-cell) interference. Notice that we made 
the choice of including the useful signal S u in h n t,u, and, as a consequence, it has to be distinguished from the commonly 
considered own-cell interference. 

With the above notations, we define the interference factor in u, as the ratio of total power received from other BS to the 
total power received from the serving BS: f u = I e xt,u/ hnt,u- The quantities f u , I ex t,u> and Iint,u are location dependent and 
can thus be defined in any location x as long as the serving BS is known. In downlink, a coefficient a, may be introduced to 
account for the lack of perfect orthogonality in the own cell. 

In this paper, we will use the signal to interference plus noise ratio (SINR) as the criteria of radio quality:7 n is the SINR 
evaluated at MS u and 7* is the SINR target for the service requested by MS u. This figure is a priori different from the SINR 
evaluated at mobile station u. However, we assume perfect power control, so SINR = 7* for all users. 

A. CDMA system 

With the introduced notations, the SINR experimented by u can thus be derivated (see e.g. [3]): 

7* = ^ (1) 

a(I intjU - S u ) + Iext,u + Noise 

From this relation, we can express S u as: 

7* 

(a + I ext u I int ,u + Noise/ h n t,u) (2) 

1 + aj* 

and the transmitted power for MS u, F bjtl = S u /gb 

using relations Ii n t,u — PbQb.u <ind fu — Iext,uf Iint,u ^s; 
Pb,u = : , 7 " * ( aP b + fuPb + Noise/ g b . u ) . (3) 
From this relation, the output power of BS b can be computed as follows: 

u 

and so, according to Eq|3] 

p 1 7^ Noise 

b = 1 — 7; ( 7TT- (5) 
1 -E« T+h7r u (« + /«) 

B. OFDMA system 

In OFDMA, the data are multiplexed over a great number of subcarriers. There is no internal interference, so we can consider 
that a(Iint, u — S u ) = 0. Since I e xt,u = J2j^b P j93- u > trie ex P ress i° n Q can t> e written, 

Pb,u9b,u 

7u = ^ 5 , Ar ■ — (6) 

E^b p j9j,u + Noise 
so we have ^ 

T« = ~ i Noise 



fu + 



Pb,v.gb,u 



Since p& << f u , typically for cell radii less than about 1 km, we can negligt this term and write 

lu = i (8) 

Ju 

Remark 

We showed that the interference factor is a key parameter to analyze wireless networks. In the aim to propose a closed form 
formula of /, we first develop a physical model of the network. 



III. Physical Analysis 



The physical network model considers that a discrete set of BS of a cellular network can be replaced by a continuum of 
BS. This approach can be compared to other approaches in physics, where discrete entities are considered as a continuum. For 
example, to analyze certain types of problems concerning the gravitation, a discrete mass M can be considered as a continuum. 
And in electrostatics, some discrete charges can be considered as a continuum, too. Hereafter, we remind these approaches 
based on the Gauss theorem, and we propose an analogy with a wireless network. We first remind the Gauss theorem, applied 
to mechanics and electrostatics, which allows to replace a discrete mass (or charge) by a continuum of mass (or charge). 
Afterwards, writing the local expression of the Gauss theorem, we remind its equivalence with the fact that the gravitational 
(electrostatic) field can be interpreted as the response of the space to a mass (charge) solicitation. We go on, observing that 
it is also the case for a BS. More precisely the power received at a given point of a wireless network can be interpreted as 
the response of this point to a BS transmitting power solicitation. And we show that a BS can be replaced by a density of 
BS, when the power received has the following form Kr~ r > with n = 2. Considering a multidimensional analysis of wireless 
networks, we show that this results can be generalized for 77 > 2. Finally we interpret that last analysis, in the case of a realistic 
network. 

A. Gauss theorem: global expression 

In a mechanical system, let's consider a mass m inducing a gravitational force at a given point of the space, located at 
the distance r from the mass. The gravitation is characterized by a function form \. The Gauss theorem shows that the 
gravitational effect of a non punctual object with a mass density p m , a volume V and a mass m can be replaced by an 
equivalent one with the same mass which is located at its gravity mass centre. For example the Earth can be replaced by 
a point in its centre where all the Earth mass would be concentrated. The gravitational effects of these two approaches are 
equivalent. In fact the Gauss theorem expresses that a field created by a non ponctual object of mass m(gravitational case) or 
charge (/(electrostatic case) at a given point of a system can be calculated by considering only the distance to the mass centre 
(or the charge centre), and not the distances to each elementary mass (or charge). The Gauss theorem allows replacing a mass 
distribution density (or charge distribution density) by one unique mass (or charge). 

1 ) Electrostatic case: The electric field E flow through a closed surface S is proportional to the sum of the charges which 
are inside this surface. 

E.dS=— (9) 
s £ o 

where Q = J J J y p c dV and p c is the density of charges in the volume V. 

2) Gravitational case: The gravitation force F flow through a closed surface is equal to the sum of the interior masses rrii 
multiplied by — AttG. 



F.dS = -4nGj2 m i ( 10 ) 



where 2~2 m i ~ f f fv PmdV and p m is the density of mass in the volume V. 

3) Application to a Wireless Network: The question is: can we replace a set of discrete transmitting sources (BS) by 
a continuum of sources (pes)? We notice that a source BS is characterized by its transmitting power P t . In other terms, 
considering the magnetic field B and the Poynting vector P = E AB characterizing the transmitting electromagnetic power, 
is it possible to write such a relation: 



where 2~2^t — J J J v PpdV and p p is the density of power in the volume V, P r is the module of the Poynting vector received 
at a given point, and P t the power transmitted by a BS? 

Since a BS is characterized by its transmitting power P t , that question means: is it possible to replace a punctual BS by a BS 
density all over a cell? We can express differently that question as follows: Is the effect of a punctual BS at a given point of 
a network the same as the effect of a BS density all over a cell (analogy with mechanics or electrostatic)? To answer that last 
question, we will express the Gauss theorem differently. Precedently, we expressed the global (or integral) expression of the 
Gauss theorem. Hereafter we analyze its local expression. 



B. Gauss theorem: local expression 

Hereafter, we show that, from the local expression of the Gauss theorem written as the Poisson equation, we can interpret 
the gravitation (and the electric) field as the Response of the Isotropic Space to a mass (charge) density solicitation [10]. 

The force F(r, t) created at each point of a system where there is a density of mass is obtained as the sum of all the 
infinitesimal contributions p(r',t) 

8F(r,t)=Gp(r',t)6 3 r'^^ (12) 

of a continuum distributed matter: 

F(r.t) = / Gp(r',t)S 3 r' I - C d 3 r'. (13) 
J ||r'-r'|| 3 

Introducing a gravitational potential ip(r,t): 

i/>(r,t) = - [ Gp(r',t)S 3 r' 1 _ d 3 r' (14) 
J \\r — r'\\ 



and when we notice that 



9rad{ ^\? = ii^iJ' (15) 
1 1 r* — r* 1 1 II r — r'\\ s 

we can write F(r,t) = —grad(4>(r,t)) , and finally: We can write the potential as a convolution: 

i>{r,t) = -Gp{r,t)*^ T (16) 

IMI 

And so grad(ip(r,t)) = —Ggrad(p(r,t) * tAt) . 

The function •jAr) is the Green function of the Laplacien A, to a constant. We have A = — 4Tr8(r) where S(r) is the Dirac 
distribution. So we have 

A(ip(r,t)) = 4nGp(r,t) * S(r) (17) 
Since the Dirac is the neutral element of the convolution algebra, we finally obtain the Poisson equation: 

A(V>(r,t))=47rGp(r,t) (18) 

These equations mean that Gravitation can be interpreted as the Response of the isotropic Space to a density mass 
solicitation. 

Concluding remark 

It is equivalent to say that the gravitation is the response of the isotropic space density mass solicitation and to write the 
equation 17 which expresses that it exists a gravitational potential ip whose Laplacien is proportional to a mass density . 

C. Equivalent of Poisson Equation for a Base Station 

We consider a BS, and the zone covered by it denoted a cell. The BS transmits the power P t . The power received at a given 
point of the network can be expressed as proportional to ^f. We can also say that each base station has an impact to a point 
of the network given by that expression. This effect can be interpreted as the Response of the environment to the BS power 
transmission solicitation. 

We thus can introduce a power potential tp p , and a density of BS pbs, an d write this response with an equivalent of the 
relation [17] for the case of a BS, as follows: 

A(<Mr)) = Kp BS *6{r) (19) 



Remark: due to the property of the Dirac function, a localized BS transmitting power can be expressed as a localized 
density of BS transmitting power over the transmission surface. 



In the case of 77 = 2, we express the analogy with the gravitation as follows. The power P r received at each point of a 
system is due to a density of transmitting BS, and is obtained as the sum of all the infinitesimal contributions (analogy with 
( fl2] i) of a continuous base station distribution. 

SP(r, t) = Kp BS {r', t)S 3 r' .. r ~ T ' s (20) 

||r — r'W* 

As a consequence, for 77 = 2, we can derive a "Base station's Poisson equation" for the power received at a given point of a 
network. In our analysis, a BS can be defined by its transmitting power P t localized at the center of the cell. It can be written 
Pt$Bs( r )- Using the Gauss theorem, we can replace a BS by a density of BS all over the cell area covered by the BS. So the 
impact of a localized BS on the network is the same as the one due to a density of BS distributed all over the cell area. As a 
consequence the localized BS can be replaced by a BS density all over the cell. 

D. Some concluding remarks 

We recalled that the gravitational (or electrostatic field) can be interpreted as the Response of the isotropic Space to a density 
mass solicitation (Poisson). Consequently the gravitational effects of a non ponctual object of mass m and density of mass p m 
are the same as the ones induced by an object with the same mass m concentrated at its mass centre (Gauss). And finally a 
mass m can be replaced by a volume with a density of mass (Gauss). Following an analogue approach with the gravitation's 
case, we showed that the power received at a point of a system coming from a BS can be interpreted as the Response of 
the isotropic Space to a BS solicitation. So we can write a "Base station's Poisson equation". Consequently the impact of the 
transmit power P t of a BS (localised at the centre of a cell) is the same as the one induced by a density of BS all over the 
cell with the same total transmitting power. And finally a transmitting BS can be replaced by a density of BS distributed all 
over the cell. 

We notice that these developments are available due to the fact that we consider a free and isotropic space (77 = 2). 

E. Wireless network analysis 

For 77 = 2, we showed that a BS can be replaced by a density of BS. However the parameter 77 is not equal to the value 2 in 
a real environment. It is generally comprised between 3 and 4. This is due to the reflections, refractions of radio propagation 
... all phenomena which transform an isotropic and free environment to a non one. 

Remark: The equation [19] can be generalized to obtain the classical definition of the gravitational potential in a D- 
dimensionnal space [10]. So ip can be written as proportional to £> 1 _ 2 in a multidimensional space where D represents 
the dimension of the space. 

That property can be interpreted as follows: the response of the system to a solicitation depends on the dimension of the space 
considered for the system. For example the gravitational force could be written as proportional to r£ , 1 _ 1 . Generally, D=3 for 
the "normal" space. In wireless network, the values of 77 should be equal to 2 in a free space. The differences observed for the 
values of 77 can be interpreted as follows. When the space is not free, the parameter 77 characterizes a pathloss parameter, for a 
multidimensional system, with D > 3, in which the space would be free and isotropic. In this case, the response of the space 
to the power solicitation of the BS takes into account the environment, as if the set of elements constituting the environment 
(trees, cars, building) would characterize other dimensions effects. In other words, since 77 = 3, the 3 dimensional-space can be 
expressed as a "four" dimensional one. This four dimensional space (D=4) is free and isotropic (in a point of view of power 
reception), and the response to the power solicitation can be written as proportional to = \ . 

Another way to express this approach consists in saying that a non-free environment in dimension 3 with 77= 3 can become 
free in dimension 4: the environment increases the space dimensions. We can apply the Gauss theorem in this 4-dimensional 
space and replace a BS (i.e. a localized power transmission) by a continuum of BS all over the surface cell. 

We notice that considering a unique value of 77 means that the environment is isotropic 

IV. A Physical Network Model 

In this section, we first present the model, derive the closed-form formula for /, and compare it through Monte-Carlo 
simulations in a hexagonal network. 

A. OCIF formula 

We showed that it is possible to replace a ponctual BS (i.e. transmitting power) by a density of BS distributed all over 
the cell. We generalize that approach to a network. The key modelling step of the model we propose consists in replacing a 
given fixed finite number of BS by an equivalent density of transmitters which are spatially distributed in the network. This 
means that the transmitting power is now considered as a continuum field all over the network. In this context, the network is 



characterised by a MS density pms and a base station density pss [7]. We assume that MS and BS are uniformly distributed 
in the network, so that pms an d Pbs are constant. As the network is homogeneous, all base stations have the same output 
power Pft. 

We focus on a given cell and consider a round shaped network around this centre cell with radius R nw . The half distance 
between two base stations is R c (see Figure [T]). 




Fig. 1, Network and cell of interest in the physical model; the distance between two BS is 2R C and the network is made of a continuum of base stations. 



For the assumed omni-directional BS network, we use a propagation model, where the path gain, g\, u , only depends on 
the distance r between the BS b and the MS u. The power, pb :U , received by a mobile at distance r u can thus be written 
Pb,u = PbKr~ v , where if is a constant and 77 > 2 is the path-loss exponent. 

Let's consider a mobile u at a distance r u from its serving BS. Each elementary surface zdzdO at a distance z from u contains 
pBS z dzd9 base stations which contribute to I ex t. u - Their contribution to the external interference is pBszdzd6P\,Kz~ ri . We 
approximate the integration surface by a ring with centre u, inner radius 2R C — r u , and outer radius R nw — r u . 



L ext.u 



27T nR nv 

2R C - 



PssPbKz n zdzdO 



= 2npBsI } K [(2R C - r u ) 2 -" - (R nw - r u ) 2 ""] . (21) 
r\ — 2 

Moreover, MS u receives internal power from b, which is at distance r u : Iint.u = PbKr~ v . So, the interference factor 

f u = Iext.u/ Iint,u can be expressed by: 

= ^PBSTl _ ru) 2- V _ {Rnw _ ru) 2-„] _ (22) 

77 — 2 

Note that /„ does not depend on the BS output power. This is due to the fact that we assumed an homogeneous network 
and so all base stations emit the same power. In our model, / only depends on the distance r to the BS and can be defined 
in each location, so that we can write / as a function of r, f(r). Thus, if the network is large, i.e. R nw is big in front of R c , 
f u can be further approximated by: 

f{r) = 2 -^f{2R c -r)^. (23) 

This closed-form formula will allow us to fastly compute performance parameters of a cellular radio network. We notice 
that formula can be applied to any wireless system (TDMA, ad-hoc.) as long as pss represents the density of interfering 
transmitting nodes. 

B. Comparison with hexagonal model network 

In this section, we aim at comparing the physical approach presented above to an hexagonal classical one. In this perspective, 



we will compare the figures obtained with Eq 22 with those obtained numerically by simulations. The simulator assumes an 
homogeneous hexagonal network made of several rings around a centre cell. Figure [2] shows an example of such a network 
with the main parameters involved in the study. 

The comparison is done numerically by computing / in each point of the cell and averaging the values at a given distance 
from the BS. Figure [3] shows the simulated interference factor as a function of the distance to the base station. Simulation 
parameters are the following: R = 1 Km, a = 0.7, 77 between 2.7 and 4, pss = (3v / 3i? 2 /2) _1 , the number of rings is 15, 



Fig. 2. Hexagonal network and main parameters of the study. 




Fig. 3. Interference factor vs. distance to the BS; comparison of the physical model with simulations on an hexagonal network with rj = 2.7, 3, 3.5, 4. 



and the number of snapshots is 1000. Eq 22 is also plotted for comparison. In all cases, the physical model matches very well 
the simulations on an hexagonal network for various figures of the path-loss exponent. It allows calculating the influence of a 
mobile, whatever its position in a cell. 

Note that the considered network size can be finite and chosen to characterize each specific local network's environment 
(urban or country, macro or micro cells). Since the physical model allows to establish a closed form formula of the interference 
factor, it allows to do analytical studies of wireless networks such as outage probability, quality of service, capacity [9]. 

V. Conclusion 

We proposed a new framework for the study of cellular networks called the physical network model. The key idea of the 
physical model is to replace the discrete base stations entities by a continuum of transmitters which are spatially distributed in 
the network. This allows us to obtain simple analytical expressions of the main characteristics of the network, the interference 
factor /. We focused on CDMA and OFDMA networks. Following an analogue approach with the gravitation's case, we 
showed that the power received at a point of a system coming from a BS can be interpreted as the response of the isotropic 
space to a BS solicitation. So we can write a "Base station's Poisson equation". Consequently the impact of the transmit power 
P t of a BS (localised at the centre of a cell) is the same as the one induced by a density of BS all over the cell with the same 
total transmitting power. Therefore a transmitting BS can be replaced by a density of BS all over the cell. Generalizing that 
approach to a whole network, we established a closed form formula of the interference factor. Comparing this network approach 
to a classical hexagonal one, we showed the physical model and the traditional hexagonal model are two simplifications of 
the reality. None is a priori better than the other but the latter is widely used, especially for dimensioning purposes. That 
is the reason why a comparison is useful. However, since the physical model allows to establish a closed form formula of 



the interference factor, it allows to do analytical studies of wireless networks such as outage probability, quality of service, 
capacity. 
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